Simulations of Coulombic Fission of Charged Inviscid Drops 
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We present boundary-integral simulations of the evolution of critically charged droplets. For 
such droplets, small ellipsoidal perturbations are unstable and eventually lead to the formation of 
a "lemon" -shaped drop with very sharp tips. For perfectly conducting drops, the tip forms a self- 
similar cone shape with a subtended angle identical to that of a Taylor cone. At the tip, quantities 
such and pressure and fluid velocity diverge in time with power-law scaling. In contrast, when 
charge transport is described by a finite conductivity, we find that small progeny drops are formed 
at the tips whose size decreases as the conductivity is increased. These small progeny drops are of 
nearly critical charge, and are precursors to the emission of a sustained fiow of liquid from the tips 
as observed in experiments of isolated charged drops. 

PACS numbers: 47.65.-d, 47.15. km, 47.55.D-, 47.11. Hj 
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An isolated droplet of liquid will naturally take the 
form of a sphere in order to minimize its surface area. 
If wc now place a net amount of electrical charge on the 
sphere, there is a pressure opposing the effects of surface 
tension due to the repulsion of mutual charges. As first 
described by Lord Rayleigh [ij , there is a critical amount 
of charge Qc that can be placed on the drop before the 
sphere will become unstable to small perturbations from 
equilibrium: 

= 87r(ecri?3)i/2 (i) 

where R is the radius of the drop, a is the surface tension, 
and e is the electrical permittivity. The simplest case, 
and the one first directly observed in experiment 0, Q , 
is that of an ellipsoidal perturbation where the droplet 
evolves into a "lemon" shape, and high-speed jets of liq- 
uid carrying a significant fraction of the total charge are 
emitted from the tips. 

The disintegration of such isolated charged drops oc- 
curs in natural settings such as thunderstorm clouds and 
bursting bubbles at the ocean surface [1], industrial ap- 
plications ranging from ink-jet printing to electrospray- 
ing la , and is especially important in mass spectrometry 
d, U Charged liquid drops were also used as an early 
model for the mechanism of nuclear fission Q. Nearly 
all previous studies have looked solely at the oscillations 
around equilibrium and the limits of stability (9l-[ll|. or 
operated under the assumption of infinite conductivity 
12|-|14||. This latter assumption is especially important; 
recent numerical simulations show that bulk conductiv- 
ity controls the fine fluid jets formed in applied electric 
fields [iBl- Indeed, the specific mode of charge conduc- 
tion will affect the dynamics of the jet [l^. The most 
popular models of charge transport suppose bulk con- 
duction 15|, |l7[ , although recent evidence suggests that 



conduction along the surface of the drop is significant 



In this letter, we simulate the initial instability and 
eventual cone-jet formation for charged inviscid drops 
with total charge Qc- We consider two separate cases: 
(1) The drop is a perfect conductor. In this regime the 
"lemon" -shaped drop forms extremely sharp tips, where 
quantities such as charge density, curvature, and veloc- 
ity diverge in finite time. The tip shape is self-similar 
and conical, and the subtended angle is exactly that of 
a Taylor cone, which results from the curvature term be- 
ing present in the asymptotic balance of forces. (2) The 
charge transport is limited by a finite electrical conduc- 
tivity. We investigate the effects of bulk and surface con- 
duction of charge; in both cases the limited conductivity 
results in a blunting of the sharp cone-tip, and the even- 
tual emission of a progeny drop from the tip. Regardless 
of the method of charge transport, the progeny drop car- 
ries an amount of charge just below the Rayleigh limit 
Qc, which depends weakly on its size. 

Numerical method — We begin with an incompressible, 
axisymmetric fluid globule immersed in an incompress- 
ible fluid of infinite extent. The axis of rotation is the 
z-axis and r is the radial coordinate. The inner fluid has 
a density pi and the exterior fluid has a density p2, where 
the density ratio A = p2/pi- With this formulation we 
can consider the effects of an exterior fluid, so that the 
behavior of both bubbles and droplets can be studied, al- 
though in this letter we will focus solely on the case where 
A < 0.001. The interface also has a uniform surface ten- 
sion cr and a varying surface charge density q with total 
charge Qc- The charge conduction along the surface can 
be perfect so that q only depends on the geometry of the 
drop, or transport properties can be specified with a bulk 
conductivity k or surface conductivity 7. The flow every- 
where is assumed to be inviscid and irrotational (smooth 
and non-turbulent). With these assumptions the veloc- 
ity of the fluid v can be described by the gradient of a 
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FIG. 1. (a) Evolution of the drop from a slightly prolate 
shape to the final "lemon" shape with pointed cone tips, (b) 
Mean curvature and charge density {q*^^) at the tip 

as a function of r for an infinitely conducting droplet in a 
vacuum. The tip region (dotted circle in (a)) seen in the 
inset shows the evolution of the droplet interface into a Taylor 
cone. The quantities and g*'*' scale as t~'^^^ and t~^^^ , 

respectively. The solid and dashed lines show power-law fits 
to the data with fixed exponents, (c) Time evolution of the 
semi-cone angle for 3 different simulations. The semi-cone 
angle always asymptotically approaches Taylor's value of « 
49.3°, regardless of initial conditions, indicating a universal 
self-similar solution. The approach of the semi-cone angle 
scales roughly as r^/^, thus the choice of abscissa. 



scalar potential "v = V0, and the problem is reduced to 
solving Laplace's equation = "with time-dependent 
boundary conditions 
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2q^-2K^ (2) 



"where Km is the mean curvature and the superscripts 
refer to the exterior (-|-) or interior (-) side of the inter- 
face. In eqn.[2]and all further discussions, quantities have 
been made dimensionless "with lengths scaled by the ini- 
tial radius R of the globule. All times in the problem are 
scaled by (R^pi/a), and any units of charge are scaled 
■with 2(ecri?^)^/^ so that critical surface charge density is 
equal to unity. For simplicity, we also assume the elec- 
trical permittivity everywhere is equal to e. 

We follo"w previous boundary-integral methods 
used to accurately compute the motion of an interface 
bet"ween t'wo inviscid fluids, "with addition of the charge 
density in eqn.[2j In our formulation there is no bulk free 
charge and the total charge on the interface is conserved. 
If the conductivity is infinite, then the charge distribu- 
tion is purely geometry-dependent and the charge "will 
al"ways arrange itself so that the electric potential is 
constant. In this case we solve the fairly simple electro- 
statics problem of finding the surface charge distribution 
on a charged axisymmetric conductor, "which involves in- 
verting an integral equation. 

Ho-wever, this is a special case. When the conductiv- 
ity is finite, "we follow the method described in reference 
[23|, "with the addition of terms accounting for Ohmic 



conduction 



Dq 
Dt 



1 



TVlip + K{n■V^P)- (3) 



"where w„ is the normal velocity of the interface, Vs is the 
surface gradient operator, and the convective derivative 
operator is D/Dt = d / dt + VnH ■ V . The dimensionless 
number Pe — {Ra/x^pY^^ is the Peclet number, "where 
X is the surface diffusivity. For all simulations, "we use 
a value of Pe=1000, so that diffusion of charge is es- 
sentially negligible. The parameters K = k{R^ p/ae^y^"^ 
and F = "/{Rp/ae'^Y^^ are a measure of the bulk and sur- 
face conductivity, respectively, and each is defined as a 
ratio of time scales, t conduction /t capillary, where the con- 
duction timescale depends on the mode of conduction, 
bulk or surface. For our simulations, we varied K and F 
independently while the other was set to zero so we could 
isolate the effects of a particular transport coefficient. 

Perfect Conductors — First we will consider the case of 
a perfectly conducting drop. For simplicity, we ignore the 
exterior fluid and choose A=0. The charge distribution 
is calculated purely based on the interfacial geometry at 
each time step. After the lemon-shape is formed (Fig. 
[Tfa)), the tip continues to sharpen indefinitely: both the 
curvature and the charge density diverge. Fig.[IJb) shows 
the mean curvature and charge density at the tip 
as a function of t — tg — t, where to is a parameter chosen 
from a power-law fit of the data, and represents the point 
of divergence. The mean curvature scales as r~^/'^, and 
the charge density as r^^/"^. Although not included in 
Fig.[IJb), we find that the normal velocity of the interface 
at the tip scales as t~^^^. These scalings indicate that 
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FIG. 2. Final drop shapes at the time of progeny drop emis- 
sion for several values of the dimensionless surface conductiv- 
ity r and bulk conductivity K. The drop is surrounded by 
dilute vapor (A=0.001). The initial evolution of the shapes is 
similar to the infinite conductivity case (Fig. [Ha)) with the 
addition of a progeny drop formed at the tips. Lower val- 
ues of r as seen in (b), or K as seen in (d), lead to more 
elongated drops and larger progeny drops because the charge 
moves more slowly. The insets are zoomed-in images of the 
tip. The specific mode of charge conduction determines the 
shape of the drop and tip prior to emission. 



all terms in eqn. [2] will be important for an asymptotic 
balance. This is in contrast to the Coulombic fission of a 
viscous, perfectly conducting droplet, where the charge 
density scales as r^"-^ and the curvature as r^° ''^ 12l |. 
so that the curvature term becomes negligible as r — >■ 
and the viscous forces balance the electrostatic forces. 
Similar exponents were also measured experimentally in 
24| . which considered the formation of a spout from an 
oil-water interface in an applied electric field. 

The power-law behavior in the dynamics suggest a self- 
similar solution for the tip region. Our results show that 
the asymptotic shape of the tip region is that of a perfect 
cone, in contrast to simulations in reference |14| with a 
much lower dynamic range in space and time. Surpris- 
ingly, the semi-cone angle is exactly equal to that of a 
Taylor cone, ar ~ 49.3°. Why is this so? Taylor ^ 
assumed a steady-state solution where the electrostatic 
pressure exactly balanced the curvature pressure on the 
surface of a cone (last 2 terms in eqn. [2]). He showed 

1 /2 

that V oc rj Pi/2{cos9) in the far field, where is 
the distance from the apex of the cone, P is a Legendre 
polynomial, and 9 = n — a is the obtuse angle to the ver- 
tical (Fig. [TJc)). For a perfect conductor ip is constant, 
so Pi/2icos9) = and 9 = n — ax- For our dynamic 



case, one realizes that the velocity potential also satisfies 
Laplace's equation, and (j) should be of the same form as 
ip. Thus capillary, electric, and Bernoulli forces all bal- 
ance in the asymptotic cone shape, and the cone angle 
must be ax- This solution is apparently universal, as 
different initial simulation conditions all converge to the 
same value for the cone angle in Fig. [TJc) . 

Finite Conductivity — Next we will consider the case 
where the charge transport is limited by a finite electri- 
cal conductivity. For these simulations we use a density 
ratio A=0.001, corresponding to a liquid drop in an am- 
bient gas atmosphere. Previous experimental studies of 
charged drops 18 1 and liquid bridges [l^ suggest that 
both bulk and surface conduction are important in ionic 
solutions due to a diffuse layer near the surface of the liq- 
uid. We have simulated both cases independently to elu- 
cidate the effects that the mode of charge transport has 
on the drop dynamics. Fig.[2]shows the final shapes of the 
drop just before progeny drop emission for two values of F 
(a)-(b), and for two values of K (c)-(d). In both cases, for 
large conductivities, the final shape resembles the perfect 
conductor, with the addition of a tiny drop emitted from 
the tip. Our simulations cannot proceed past the emis- 
sion of the first drop, after which a jet /stream of droplets 
have been visualized in experiments [3, [3] . Lower con- 
ductivities result in large progeny drops and highly elon- 
gated drop shapes, which is due to the reduced rate of 
charge transport to the pointed tips. As the tip sharp- 
ens, the charge must move to the tip region to remain in 
equilibrium. If this cannot occur quickly enough, then 
the surface tension forces will act to breakup elongated 
structures through the Rayleigh-Plateau instability. At 
the final moments of progeny drop pinch-off, the charge 
density remains finite and the dynamics are governed by 
the known universal solution for inviscid capillary pinch- 
off 

Fig. El^a) shows the dependence of the progeny drop 
radius Rp and charge Qp on the conductivity parameters 
F and K. For bulk conduction, Rp oc K^"^/^ and Qp oc 
K'"^. The progeny drop size can be understood as a bal- 
ance of the conduction time scale tconducUon = e/^ with 
the local capillary time scale tcapUiary = {RpP/'^Y^^ ■ 
If we assume that the progeny drop has of order the 
Rayleigh limit of charge (eqn. [T]), then we obtain the 
K^^ scaling for Qp as well. This argument is identical 
for the case of surface conduction where Rp oc F"^ and 
Qp (X r~'^, except the conduction time scale is now de- 
fined by the surface conductivity tconducUon = ^RpH- 
The observed scalings in the simulations are in excellent 
agreement with these predictions. Fig. ^[h) shows the 
ratio of Qp to the Rayleigh limit for the progeny drop 
87r(e(T i?p)^/^. Regardless of the charge conduction mech- 
anism, this ratio only depends on the progeny drop ra- 
dius Rp, and slowly approaches unity for small drops. 
This suggests that drops emitted from the tips during a 
Coulombic fission process are only marginally stable, and 
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FIG. 3. (a) Radius of the progeny drop Rp (closed symbols) 
and the total charge on the progeny drop Qp (open symbols) 
at the moment of pinch-off. Red circles refer to simulations 
using surface conduction and blue squares using bulk conduc- 
tion. The solid and dot-dashed lines show a power-law fit 
of and r^^, respectively. The dashed and dotted lines 
show power-law fits of and K~^^^. These scaling laws 
are expected from equating the conduction and capillary time 
scales for surface and bulk conduction. As the conductivity 
is increased, the size of the progeny drops become smaller 
and the total charge they carry decreases, (b) Ratio of the 
total charge on the progeny drop to the critical amount of 
charge necessary for the instability (eqn.[T]). This is the frac- 
tion of the Rayleigh limit, so that a value of I means the 
progeny drop will certainly be unstable and undergo further 
fission. Regardless of the conduction mechanism, this fraction 
depends only on the drop size. 



could be subject to further breakup given a sufficient per- 
turbation. 

Conclusion — For the first time, we provide a quanti- 
tative picture which shows how charge conduction con- 
trols the shape and eventual emission of drops during the 
Coulombic fission of inviscid drops. For perfectly con- 
ducting drops, a self-similar cone shape is formed with a 
cone angle of « 49.3°, identical to Taylor's steady-state 
value. Although our simulations are dynamic and the 



pressure and velocity are diverging at the tip of the cone, 
this cone angle can be understood solely by the presence 
of the mean curvature in the asymptotic balance of forces. 
When the charge transport is limited by a finite surface or 
bulk conductivity, a small progeny drop is emitted from 
the tip of the cone; a precursor to the tip-jetting observed 
in experiments. The drop size and charge are determined 
by a ratio of the conduction time scale to the capillary 
time scale, and the progeny drop is nearly unstable and 
may undergo further Coulomb fission. 
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